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ABSTRACT 

Polarization has proved an invaluable tool for probing magnetic fields in relativistic jets. Maps of the intrinsic 
polarization vectors have provided the best evidence to date for uniform, toroidally dominated magnetic fields 
within jets. More recently, maps of the rotation measure (RM) in jets have for the first time probed the field 
geometry of the cool, moderately relativistic surrounding material. In most cases, clear signatures of toroidal 
magnetic field are detected, corresponding to gradients in RM profiles transverse to the jet. However, in many 
objects these profiles also display marked asymmetries which are difficult to explain in simple helical jet models. 
Furthermore, in some cases the RM profiles are strongly frequency and/or time dependent. Here we show that these 
features may be naturally accounted for by including relativistic helical motion in the jet model. In particular, we 
are able to reproduce bent RM profiles observed in a variety of jets, frequency dependent RM profile morphologies 
and even the time dependence of the RM profiles of knots in 3C 273. Finally, we predict that some sources may 
show reversals in their RM profiles at sufficiently high frequencies, depending upon the the ratio of the components 
of jet sheath velocity transverse and parallel to the jet. Thus, multi-frequency RM maps promise a novel way in 
which to probe the velocity structure of relativistic outflows. 

Subject headings: galaxies: jets - magnetic fields - polarization - radiative transfer - radio continuum: general - 
techniques: polarimetric 



1. INTRODUCTION 

Polarization observations of relativistic jets have been 
, instrumental in identif ying the presence of l arge-scale or- 

■ dered m agnetic fields ("Tavlor & Perlev* '199 2t iTavlod [1998 : 

■ [ Lister et al. 1998; Taylor 2000; Zavala & Tavlod 2001 , 

■ '2002', '2003, 2004; Asada et al.1 l2008alibt iKharb et all [2009 ; 
£)' Sullivan & Gabuzda 200^ For nearly all active galactic 
nuclei (AGN) jets, Faraday rotation dominates the polarization 

, signal, scrambling the polarization angles of the intrinsic 
emission on angular scales small compared to the jet width. 
This can be ameliorated by determining the rotation measure 
(RM), via multi-wavelength polarization observations, and 
then removing the associated distortion using the standard 
relation, 

*mt=*obs-RMA2, (1) 

where ^'im.obs is the intrinsic/observed polarization angle. Typ- 
ical RM's for AGN jets and cores range from lOradm"^ to 
lO-'radm"-, though it can reach lO^radm"^ in some extreme 
cases. Thus, at IGHz frequencies, the intervening Faraday 
rotation results in polarization angle rotations of roughly 6°- 
600°. Nevertheless, once the Faraday rotation is solved for, 
via multi-wavelength polarization observations, and removed, 
most jets show remarkably uniform polarization maps, gener- 
ally aligned with the jet axis. Both the polarization fraction and 
the orientation has been used to argue that in the observer frame 
jets are overwhelmingly toroidally dominated (iLyutikov et al.l 
6005). 

Observations of the RM itself provide information about 
the intervening magnetized plasma. There are multiple po- 
tential sources for the RM's including the Galactic interstel- 
lar medium, the intracluster medium near the AGN and cooler, 
less relativistic material surrounding the jet itself. The Galac- 
tic RM itself can be on the order of 10^-10^ radm"^ and varies 
considerably on large scales, presumably as a consequence of 
the structure of the Galactic magnetic field dEkers et al.ll 19691; 



iMen et alj l2008t iNoutsos et all l2008h . The magnetic field 
strength and geometry within the intracluster medium is poorly 
constrained, though few sources lie sufficiently near th e centers 
of co oling core clusters for this to dominate (Carilli & Taylor! 
I2OO2I) . In contrast, Faraday rotation intrinsic to the jet environ- 
ment itself is expected to produce large RM's. 

The morphology of the Faraday rotating medium surround- 
ing AGN jets is still poorly constrained. While the lack of 
Faraday depolarization implies that the medium must be in the 
foreground and itself produce negligible emissions, presently 
it is not possible to distinguish between nearby unrelated ion- 
ized gas clouds and rotation within the near side of a jet sheath, 
presumably containing an ordered helical magnetic field. Free- 
free absorption from pc-scale foreground clouds has been de- 
tected in Cen A and NGC 1275, implying that these clouds 
exist in at least some sources (Jones et al. 1996; Wal ker et al.l 
I2OOOI) . In these one would expect roughly random RM gra- 
dients, associated with the random orientations of the cloud 
magnetic fields. However, in a handful of cases thus far or- 
dered RM gradients transverse to the jet axis have been ob- 
served, providing tantalizing evidence for the presence of he- 
lically m agnetized sheaths surrounding AGN jets (Asada et al.l 
l2008allbHaiarb et alJl2009l: lO' Sullivan & Gabuzda 2009) . It is 
this latter possibiUty that we restrict ourselves here. 

For static, axisymmetric jet sheaths, the transverse RM pro- 
files are roughly linear and symmetric about the jet axis, up 
to a uniform offset depending upon the mag netic pitch angle . 
While for m any sources this is sufficient (Asada et al.ll2008bl : 
IKharb et all 120091: lO'SuUivan & Gabuzdal l2009l) . in a num- 
ber of cases the RM profiles display significant asymmetries 
( Asada et al. 2008b; O' Sullivan & Gabuzda 2009). Further- 
more, these asymmetries can be both time and frequency de- 
pendent. While relativistic motion along the jet can reduce the 
RM in the observer frame considerably , it is unable to produce 
these asymmetric profiles (lL6pezll2006l) . Here we consider the 
effects of relativistic helical motion. This is expected theoret- 
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ically in both magnetohydrodynamic and force-free jet mod- 
els, as well as motivated by jet simulations, though these typi- 
cally cannot probe the necessary dynamic range to make direct 
com parisons ( Vlahakis & Konigl 2004; De Villiers et al. 2005; 
McKinnev 2006; Komissarov et al. 2007; Tchekhovskov et al. 
2008; McKinney & Blandford 2009). Within the context of a 
simple helical jet model, it is possible to reproduce both, the 
observed asymmetries and the frequency & temporal behavior 
of jet RM profiles, without appealing to asymmetric structures 
or special viewing angles. 

2. FARADAY ROTATION IN RELATIVISTIC BULK FLOWS 

Faraday rotation is a result of the different phase veloci- 
ties of the two, nearly circularly polarized', electromagnetic 
eigenmodes of magnetized plasmas. The subsequent accrual 
of a phase difference results in a local rotation of the polariza- 
tion plane. Critical to this interpretation is the assumption that 
the plasma eigenmodes propagate nearly adiabatically. While 
changes in the underlying plasma properties (density and mag- 
netic field) can lead to violations of this condition, for simplic- 
ity we will only consider cases in which the relevant plasma pa- 
rameters change sufficiently slowly that these may be ignored 
dlnBlin 2003; Broderick & Blandford 2009). Furthermore, we 
will assume that the plasma smoothly joins the sub-relativistic 
material surrounding the jet, and thus will not consider the rel- 
ativistic aberration of the resultant polarization. In this case, 
locally, the polarization angle, 5*, changes by an amount iden- 
tical to the phase difference accrued, That is, in the plasma 
rest frame 

167r^e^ ~ ~ 
dvE' = d$=^=-^nB-d£, (2) 
nf^c V 

where n and B are the proper electron density and rest-frame 
magnetic field, v is the rest-frame frequency and dl is the rest- 
frame distance element. Since the total accrued phase is a 
Lorentz scalar, and we assume the plasma smoothly joins a 
non-relativistic flow, we obtain 

A*=l^/^B.d£. (3) 

Presently, this integral is performed in the comoving frame 
of the plasma. However, we may use the scalar nature of ^ 
to rewrite it in terms of the observed frequency, v, and an 
integral in the observers frame, Al. To do this, let us begin 
by defining the plasma four- velocity^, = 7(1, /3). In either 
the force-free or magnetohydrodynamic prescriptions, the elec- 
tric field vanishes in the plasma frame, and thus F'^^m^ = 0, 
where F^^'^ is the electromagnetic field tensor This implies that 
ptiv _ ^P-'^°'f^ii^iyp/2, where b^^ = (b' ,h) is the four-vector cor- 
responding to the magnetic field in the plasma rest frame and 
orthogonal to m^. Due to the latter condition, the components 
of are not independent, and in particular b' = /3 -b. Finally, 

' This is generally true only for cold, ionic plasmas. In pair plasmas the 
eigenmodes are lineai; while for plasmas with relativistic temperatures or sub- 
stantial non-thermal components, the eigenmodes are typically significantly 
elliptical. However, the detection of RM gradients across AGN jets provides 
strong evidence for the presence of a significant, if not dominant, ionic com- 
ponent in the plasma surrounding the jet. For this reason, we have ignored the 
possible presence of a pair plasma in our analysis. Should the jet sheath be 
comprised of an admixture o f the two, n simply represents the excess electron 
density (Hall & Shukla'2005). That the plasma is sufficiently cold is generally 
a good assumption far from internal shocks, and observationally supported by 
the lack of significant Faraday depolarization, and hence emission from within 
the Faraday rotating sheath. 

^ Henceforth we set c = 1 unless otherwise noted. 



let us define the electromagnetic wave-vector, = 27rj/(l,k), 
which is related to the electromagnetic wave trajectory, x^{rj), 
by djc^ /drj = k^^ . This implies that the line element, as measured 
in the plasma frame is d^^ = dx^ + Uiydx^u^^ = (k^ + Ui^k'^u^)dri , 
and thus, d£ = -M^fe^'d?]. Similarly, in the lab frame, d£ = k'drj. 
Therefore, 

^ = -^ = 7(l-/3-k)- 

B-di = /7^d/^= (k-/3) -bdf. 

Inserting these into the expression for rotation of the polariza- 
tion gives, 

, , 167r3e^ 1 f (k-/3)-b 

The associated rotation measure is then simply 

aA* f (k-/3)-b , 
RM= —-^ = 0.812 / yd£radm-2, (6) 

J 72 (1-/3 -k) 

where n, b and t are measured in cm"^, iiG and pc, respectively. 

When /3 = 0, this rotation measure reduces to the standard 
expression, with n and b now corresponding the plasma den- 
sity and magnetic field strength in the observer frame. How- 
ever, when /5 y 0, some care must be taken in interpreting these 
quantities, as both are defined in the plasma rest frame. In par- 
ticular, the observer frame magnetic field, B, is related to b by 

B = 7(l-^/3)-b ^ b= - (l + 7^;3/3) -B, (7) 

7 

which may be found by inspecting directly. For large 7, 
B will appear dominated by the components transverse to the 
plasma motion. 

The relativistic motion of the plasma has two important con- 
sequences. The first is to change the magnitude of the magnetic 
field along the line of sight due to relativistic aberration, pro- 
ducing the velocity dependent term in the numerator of Equa- 
tion (|6]l. The second is simply to Doppler-shift the relevant 
plasma frequency, signified by the denominator in Equation 
(|6]l, resulting in a correspondingly higher or lower RM as mea- 
sured in the observer frame. This second effect can produce 
profound variations in the RM when /3 passes nearly parallel 
to the line of sight. The ratio of the RM of approaching and re- 
ceding plasma flows is roughly 47^, which for even moderate 
Lorentz factors can be considerable. 

3. ROTATION MEASURE PROFILES OF JETS 

Having an expression for the relativistic rotation measure in 
hand, we now must specify the properties of the Faraday ro- 
tating medium. For illustrative purposes we adopt an extraor- 
dinarily simplified model: a cylindrically symmetric plasma 
sheath surrounding a jet core, appropriate at large distance from 
the jet launching region. While the cylindrical approximation is 
poor for viewing angles within the jet opening angle, for most 
cases of interest it is sufficient, allowing us to remove a degree 
of freedom from the modeling of the Faraday screen. The den- 
sity peaks at some distance from the jet axis, falling exponen- 
tially beyond (our results are qualitatively independent of the 
radial density profile). We assume helical geometries for both 
the magnetic field and the plasma four-velocity. Explicitly, in 
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Fig. 1 . — Rotation measure profiles of a static jet (u'.u^ = 0) as viewed at 
= 60°, for magnetic pitcli angle, ai,, ranging from 10° (red) to 90° (blue) 
in steps of 10°. Top: The RM profile normalized such that the variation in the 
rotation measure is unity. Bottom: The RM profiles in units of 81.2radm"^, 
con'esponding to the RM generated by Faraday screen with electron density 
Icm"'', magnetic field 1 /tG and line-of-sight depth Ipc. In both panels the 
grayed regions show areas outside the jet core in projection, and thus presum- 
ably in which the RM cannot be measured. 



cylindrical coordinates, 



n = no 



I'Jet 



= (3 cos a,, , M* = , u'^ 



bosinai, 
' R 

/3sina„ 



(8) 



where R & 4> are the cylindrical radius and polar angle, re- 
spectively, Wjet is the characteristic width of the jet sheath, no 
and bo are the density and field strength scaling (and are ul- 
timately unimportant to the question of the RM asymmetry), 
P is the velocity and a,, and ab are the velocity and magnetic 
field pitch angles, respectively. Finally, we define the angle be- 
tween the line of sight and the jet axis to be 8, i.e. k- = cos 8, 
k'^ = sin 8 cos and k'>' = -sin 8 sin Thus, 8 = 0° and 
8 = 90° correspond to lines of sight along and orthogonal to 
the jet axis, respectively. For a variety of values for these pa- 
rameters, we compute the RM as defined in Equation (|6]l as 
a function of transverse position, thereby constructing profiles 
across the jet. 

The spatial variations in the density and orientation of the 
magnetic field necessarily lead to variation in the RM across 
the jet. In the limit of a static jet, where the plasma is at rest 
in the observer frame, the RM profile is shown in Figure [T] for 
a variety of magnetic pitch angles, as viewed from 8 = 60°. 
For purely toroidal magnetic fields (ab = 90°, blue), the RM 
profile is symmetric and approximately linear, only departing 
from a line at the sheath boundaries due to the precipitous de- 
crease in the electron density. In contrast, for magnetic fields 
nearly parallel to the jet (ah = 0°, red) the background density 



Fig. 2. — Rotation measure profiles of an ultra-relativistic, purely outflowing, 
jet, M- = 10 and o„ = 0°, as viewed at © = 60°, for magnetic pitch angle, aj, 
ranging from 10° (red) to 90° (blue) in steps of 10°. Top and bottom panels 
are defined as described in Figure[T] 

variation is clearly imprinted in the RM profile. Due solely 
to geometry, the degree to which this occurs depends upon 8, 
with the largest deviations at small 8 and ab- However, the 
absolute variation in the RM across the jet in this case is con- 
siderably reduced, and in the limit of ab = 0° vanishes com- 
pletely. More importantly, polarization maps of jets typically 
imply Rb'I'/b' > 1 (Lyutikov et al. 2005). 

Jet sheaths with relativistic bulk motion along the jet axis, 
but without any helical motion, are similar to their static coun- 
terparts. In this case the primary consequence of the motion for 
the RM profiles, apart from the drastic reduction in the net RM, 
is relativistic aberration, which effectively results in observers 
at oblique angles viewing the jet from behind. This is seen 
explicitly for an extreme case in Figure |2l in which the RM 
profiles are nearly identical to those for the static jet though 
rotated 180°. Most importantly, the sense of the RM gradient 
does not reverse as a consequence of the bulk motion^. The 
RM profiles do not significantly differ from the case shown for 
7 > 3, and for smaller velocities become even more degener- 
ate in ai,. Therefore, generally, within the context of axially 
symmetric jet models, magnetic fields with moderate pitch an- 
gles are incapable of producing strong asymmetric features in 
the RM profiles, though ab and jet Lorentz factor do play a 
substantial role in determining the absolute RM. 

In stark contrast, even moderately relativistic helical motion 
easily produces dramatic asymmetric features in RM profiles, 
independent of the magnetic field pitch angle, and in some 
cases reversing the sense of the RM gradient. As seen in Fig- 
ure|3] this occurs both for trans-relativistic and ultra-relativistic 
jets, and for velocity pitch angles as low as 30°. Generic of 
relativistic helical motion is the noticeable bow in the RM 
profiles. This feature is only weakly dependent upon and 
present for a„ > 30° for u~ > 1. Unlike similar features for 
non-rotating jets, this occurs for a wide range of absolute rota- 

^ For bulk motion along the jet to reverse the RM gradient would require 
© > 7"' while ai, < 7"' . While the former is almost certainly true, the latter 
would require magnetic field geometries that are nearly parallel to the jet axis. 
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Fig. 3. — Rotation measure profiles of jets with helical velocity fields for u~ = 0.5 {top left), u~ = I (top right), ir = 3 {bottom left) and ir = 10 {bottom right). In all 
cases, ai, = 90°, = 60° and a„ ranges from 0° (red) to 80° (blue) in steps of 10° for the solid lines. Dotted fines show the RM profiles for a„ distributed in steps 
of 0.1/7 about 0, highlighting the transition due to relativistic aberration. Top and bottom panels are defined as described in Figure[T] We choose to fix if, and not 
7, since u~ is directly measured for the jet core by observations of superluminal motion. 



tion measures. 

When the revolving plasma in the jet sheath is approaching 
almost directly, i.e., /3 is within an angle 7"' of k, the direc- 
tion in which relativistic aberration rotates the magnetic field 
rapidly changes (see the dotted lines in Figure O. As a result, 
the RM evolves rapidly with a„. At the same time, the Doppler 
shift reaches its maximum, producing dramatic enhancements 
in the absolute RM, and for a„ very close to 8, complex RM 
profiles. However, for all but the slowest jets, this is restricted 
to such a small regime in viewing angle that it is relatively un- 
likely to be observed in practice. A much more important con- 
sequence is the reversal of the sign of the RM profile gradient 
at a,, = 8. This is seen most clearly in the bottom two panels 
of Figure [3] in which the RM profiles associated with velocity 
fields with high pitch angles (blue) are nearly vertical reflec- 
tions of those with low pitch angles (red). Thus, not only the 



symmetry of the RM profile is dependent upon m"^, but also the 
direction of the RM gradient itself. 

4. DISCUSSION 

The presence of relativistic helical motion in the Faraday ro- 
tating sheaths of AGN jets has a variety of specific and robust 
observational consequences, some of which have already been 
observed. These include the presence of significant, frequency 
dependent asymmetric features in the RM profiles, and even 
reversals in the RM gradient. These are illustrated for a sim- 
ple jet sheath model in Figure |4] In this case, u~ is constant 
along the jet while u'^, and thus a„, decreases with height. Low 
frequencies probe regions far from the AGN core (see, e.g., 
Blandford & Konigl 1979), where the transverse velocities are 
low and RM profile is approximately linear Conversely, in- 
termediate frequencies probe regions closer to the AGN core. 
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Fig. 4. — Rotation measure profiles of an evolving jet. The dotted lines show 
the envelope of the Faraday rotating sheath while the helical solid line is the 
path taken by a particle as it propagates within the sheath. Throughout the jet 
ir = 3. Below the red bar, a„ > (again taken to be 60°), while above a„ < 0. 
Positions along the jet where the RM profiles are shown at the right are marked 
by green bars. The RM profiles are colored as in Figure[3] Since the variation 
of the RM is dependent upon the magnetic field strength and plasma density 
in the jet, only the ratio RM/ARM is shown. In each plot of the profile, the 
components of the four- velocity parallel and transverse to the jet axis are listed. 

where the transverse velocities are higher, which for even mod- 
estly relativistic motions results in the characteristic bent RM 
profile. Sufficiently high frequencies may probe regions for 
which Qfj, < Q, if they exist, and thus exhibit reversed RM pro- 
files; 

lO'Sullivan&Gabuzdal (120091) present RM profiles for six 
blazars, as measured in a variety of frequency bands. In at least 
two cases the RM profiles display the bent shape characteris- 
tic of helical motion: 1156H-295 and 1418H-546. In 1418H-546 
this is seen clearly in a RM profile constructed from observa- 
tions between 4. 6 GHz and 8.9 GHz. At higher frequencies, the 
RM profile shows evidence of being reversed, though the ab- 
solute RM is reduced by nearly two orders of magnitude, again 
consistent with the presence of relativistic motion. In contrast, 
at 4.6-8. 9GHz, 1 156H-295 exhibits a nearly linear RM profile, 
while at 7.9-15.4GHz, the RM profile becomes dramatically 
bent. At higher frequencies (12.9^3 GHz), it becomes roughly 
linear again, which may be evidence of reduced at small 
radii, a particularly fortuitous view of the transition to a re- 
versed profile, or failure to fully resolve the RM gradient. In 
addition 0954-1-658 and 2200-1-420 present dramatic changes in 
the absolute RM, including its overall sign, with changing ob- 
serving frequency interval and radial distance from the AGN 
core. These are consistent with the observation of regions tran- 
sitioning between ultra-relativistic and moderately relativistic 
helical motion close to the AGN core. 



RM observations of knots in 3C 273 bv lAsada etall (l2008ah 
provide extraordinary examples of the bent profiles character- 
istic of helical motion (cf. Figure 4 of that paper). More inter- 
estingly, the RM profiles have evolved from 1995 to 2002, with 
the absolute RM increasing by roughly 40% and the location of 
the bend in the profile moving towards the jet center. This be- 
havior is consistent with moderate slowing of the sheath during 
this period, i.e., a decrease in the sheath velocities of roughly 
20% as they propagated outwards. 

Unfortunately, most jet RM observations only marginally re- 
solve the transverse jet structure. As a consequence it may be 
premature to distinguish between Fara day rotation in ne arby, 
though unrelated, gas clouds (see, e.g.. lLaing et alJl2008h and 
the ordered, helical sheath models discussed here. We note that 
3C 273 is an exception, is clearly resolved and provides a strik- 
ing agreement with th e profiles we present. However, space- 
based VLBI (VSOP-2 iHagiwara et al.ll2009h and millimeter- 
VLBI (e.g., iBroderick & Loebl 120091; fPoeleman et all l2009l) 
both promise order of magnitude improvements in resolution, 
and have the potential to resolve the nature of AGN jet Faraday 
screens conclusively. 
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